Particle dynamics in sheared granular matter 
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The particle dynamics and shear forces of granular matter 
in a Couette geometry are determined experimentally. The 
normalized tangential velocity V{y) declines strongly with 
distance y from the moving wall, independent of the shear 
rate and of the shear dynamics. Local RMS velocity fluctua- 
tions 5V{y) scale with the local velocity gradient to the power 
0.4 ±0.05. These results agree with a locally Newtonian, con- 
tinuum model, where the granular medium is assumed to be- 
have as a liquid with a local temperature SV{y)^ and density 
dependent viscosity. 



An important property of granular matter is partial 
fluidization in response to shear stresses 0]. Station- 
ary granular matter can sustain normal loads and shear 
stresses, but if a threshold shear stress is exceeded, part of 
the material starts to flow with properties that appear to 
differ from those of a Newtonian fluid. Unlike ordinary 
fluids, granular materials do not exhibit intrinsic ther- 
mal motion. Instead, the granular 'temperature', gener- 
ally defined as the square of RMS velocity fluctuations 
SV^, is created by the flow itself. As a result, the mean 
flow and RMS fluctuations are related. This fundamen- 
tal connection has been investigated experimentally 
but remains poorly understood. 

The flow of sheared granular materials has been in- 
vestigated in the steady state in several experiments by 
shearing material in a Couette cell between a station- 
ary outer cylinder and a rotating inner cylinder 
All experiments indicate that the mean particle velocity 
parallel to the shear direction V{y) decreases faster than 
linearly away from the inner cylinder. 

The velocity profile in three dimensions was deter- 
mined by Mueth et al. Measurements were carried 
out both in the interior of the material using X-ray and 
NMR techniques, and on the bottom surface of the Cou- 
ette cell by optical imaging. These measurements showed 
that the velocity profile on the bottom surface and in the 
interior are the same. V{y) varies from nearly Gaussian 
for kidney shaped particles to roughly exponential for 
spherical particles. Layering of particles is suggested as 
the cause of this variation. 

In previous studies in a planar geometry we have 
found that most of the fiow is confined to 5 — 6 layers 
of particles close to the shear plane. The mean parti- 



cle velocities during brief slips of the shearing plate de- 
crease roughly exponentially with distance away from the 
moving plate, consistent with the findings in the Couette 
geometry. 

The aim of the present experiment is to understand the 
relationship between mean velocities, RMS fluctuations 
and shear forces of a steady state shear flow. We have 
developed a locally Newtonian, continuum model that 
describes the granular medium as a liquid with nonuni- 
form temperature and density dependent viscosity. The 
interplay between mean flow and RMS velocity fluctua- 
tions can be understood quantitatively in this context, 
as we demonstrate. 

In order to accomplish fluidization independent of 
shear, we apply an upward airflow at a variable rate 
through granular matter sheared in a Couette geometry. 
We measure both the mean particle velocities V{y) and 
the velocity fluctuations 6V{y) on the upper surface of 
the granular material. These should approximate parti- 
cle motion in the interior based on the previous 3D mea- 
surements described above Measurements of shear 
forces in air fluidized granular matter, including a dis- 
cussion of previous related studies of shear forces Q , will 
be presented elsewhere 

In the experimental apparatus the granular material 
(0.75 mm diameter black glass beads) is confined to a 
15 mm gap between a stationary outer cylinder and a 
rotating inner cylinder (r = 51 mm), as shown in Fig 
The inner cylinder is hollow to reduce its inertia and is 
coated with a monolayer of randomly packed glass beads. 
The outer cylinder is made of glass and is coated with 
a monolayer of randomly packed glass beads up to the 
height of the top surface, which allows observation of the 
top layer of grains through a mirror as shown in Fig |^. 
To shear the material, the inner cylinder is rotated at a 
variable rate of 0.001 — 1 Hz. The lower 38 mm of the in- 
ner cylinder is stationary in order to minimize boundary 
layer effects. 

The inner cylinder is connected to a microstepping 
motor via a flexible tempered steel spring. The spring 
bending is proportional to the applied shear force. We 
measure the spring displacement with a capacitive dis- 
placement sensor that is rigidly connected to the motor 
shaft. This spring configuration allows us to measure 
instantaneous shear forces with excellent dynamic range 
and precision, and permits both stick-slip dynamics and 
continuous motion of the inner cylinder, depending on 
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parameters. The trajectories of roughly 100 individual 
particles in the surface layer are determined with a fast 
CCD camera at 30 — 1000 frames/sec. Particle motion is 
extracted from sequences of 2000 images with a spatial 
resolution of < 0.1 pixels. From the particle tracks we 
determine average particle velocities V{y) and RMS ve- 
locity fluctuations SV{y) as a function of distance y from 
the rotating inner cylinder. 




FIG. 1. Experimental setup: The granular material (be- 
tween two concentric cylinders) is fluidized by an upward air- 
flow and sheared by rotation of the inner cylinder, which is 
connected to the motor through a flexible spring (S). Shear 
forces are determined from the spring displacement. Particle 
motions in the top layer are measured through the glass outer 
cylinder with a fast CCD camera. 

The shear dynamics without airflow are found to be 
very similar to the dynamics of a plate sliding across 
a granular layer At low shear rates the motion of 
the inner cylinder is intermittent with short, rapid slips, 
and long periods of sticking. At sufficiently high shear 
rates (or with a stiff spring), steady motion of the inner 
cylinder is observed. 

Airflow has a strong effect on the shear forces and shear 
dynamics. Without airflow, stick-slip motion is observed 
with a mean shear force of 95 N/m^. At small airflow, 
the shear force is reduced, but stick-slip motion persists. 
At high airflow rates (with a mean air speed of 1.6 m/s 
through the granular material), shear forces are reduced 
by a factor of 4 to 20 N/m^ and steady sliding motion 
is observed. The experiments described below were car- 
ried out at large airflow just below the threshold of the 
'bubbling' instability. For 0.75 mm diameter glass parti- 
cles no motion on the upper surface is observed at this 
airflow in the absence of shear. The mean shear stress 
does not change with shear rate over more than two or- 
ders of magnitude in the shear rate. The shear forces will 
be discussed in detail elsewhere Here we note that 
airflow strongly reduces the shear strength of the mate- 
rial and changes the dynamics from stick-slip motion to 



steady sliding. 

The average velocity of particles at the surface of the 
shear cell is shown in Fig. ^ as a function of distance y 
from the rotating inner cylinder, for inner wall velocities 
U ranging from 0.004 to 0.4 Hz. As found previously, the 
velocity profile decays strongly to zero over a few particle 
diameters d. The normalized velocity profile is indepen- 
dent of the imposed shear rate over at least two orders 
of magnitude in shear rate. We also find that without 
airflow (at U = 0.01 Hz, solid triangles), the inner cylin- 
der moves with stick-slip dynamics rather than steady 
sliding, but the velocity profile remains unchanged. 
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FIG. 2. Mean particle velocity (normalized by the shear 
rate) as a function of distance from the inner cylinder (in 
particle diameters). The respective shear rates U (in Hz) are 
: 0.004 (hexagons), 0.04 (squares), 0.01 (open triangles), 0.4 
(crosses). The solid triangle is the velocity profile at (7 = 0.01 
Hz without airflow. The normalized velocity profile is inde- 
pendent of shear rate and shear dynamics (intermittent or 
steady motion). The dashed line is the solution of eqs. (2) 
and (3), with 5 = 4.7d, j/„ — 2.8d, and a = .4 ( see text for 
details) . 

Figure ^ shows the RMS velocity fluctuations, which 
have not been previously measured in a 3D system. The 
fluctuations are slightly larger in the direction parallel 
to the mean flow than perpendicular to it. Since paral- 
lel fluctuations would also include the effect of fluctua- 
tions in the mean flow, we show only the perpendicular 
fluctuations. The velocity fluctuations decrease roughly 
exponentially far from the inner cylinder. However, the 
fluctuations decrease more slowly with y than does the 
average velocity. 

The RMS fluctuation is the key quantity in a flow- 
ing granular material. As already proposed by Reynolds 
]lO| , the system has to dilate in order to allow flow, which 
implies particle motion transverse to the flow direction. 
This transverse motion vanishes on average and, there- 
fore, its RMS fluctuations contain the physically relevant 
information characterizing the flow properties. 
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FIG. 3. RMS velocity fluctuations perpendicular to the 
shear direction. Fluctuations decrease roughly exponentially 
far from the inner cylinder, but more slowly than the mean 
flow. The rms fluctuations are rescaled (shifted vertically) 
such that all epxerimental points flt on the same master curve. 
The dashed line is the theoretical result (see text), with a de- 
cay length 5 — 4.7d and a boundary position = 2.8d. 

In a phenomenological description, it is tempting to 
consider these RMS velocity fluctuations as an effective 
granular temperature and consider a local hydrodynamic 
model as a starting point, as has already been proposed 
in the literature (see e.g. ||ll|). Within such a picture the 
granular temperature obeys a transport equation charac- 
terizing the balance between "heat flow" and dissipation 
due to the inelastic collisions : 



dyXdyT - eT = 0. 



(1) 



(Note that a local heating term has been omitted, which 
can be shown to be valid far from the moving boundary.) 
This local equation introduces two transport cocfhcients, 
the thermal conductivity A and the energy loss rate e. 
Both A and e are proportional to the collision frequency. 
In the high density limit we consider, p ~ with 
the density at random close packing, A and e are propor- 
tional to P/T^^^, where P is the isotropic pressure [pd]| . 
Stress conservation in the y direction implies that P is 
constant over the cell. From eq. (Q), one can obtain the 
temperature profile as 



_/2 Cosh{{H -y)/6) 
coshiH/S) 



(2) 



where S, defined as (5^ = 2A/e, is proportional to the bead 
diameter d, and is independent of position and tempera- 
ture. A vanishing heat flux has been assumed at the outer 
stationary boundary in agreement with additional exper- 
imental measurements for wider shear regions, while the 
temperature at the moving boundary, Tq, is introduced 



as a boundary condition at the wall. Tq is created by the 
flow in the flrst layers close to the wall, so that eq. is 
assumed to apply only for distance y larger than a cutoff 
distance y^ from the wall : y^ quantifies the width of 
the boundary layer in which temperature is created by 
the motion of the inner cylinder. The fit to experimen- 
tal data yields y^ = 2.8d. The effective cell width H is 
thus Hq — y-w, where Hq is the wall to wall distance. As 
shown in Fig. ||, this hydrodynamic model shows good 
agreement with the experimental data. 

Within the framework of this phenomenological model, 
the velocity profile can be constructed from the con- 
stant shear-rate condition, a^y = vi — const, with 
7 — dVx{y)/dy the local shear rate and rj the viscos- 
ity. The viscosity is expected to scale with the col- 
lision frequency, rj = rjoP/iPcd^T^^'^) (where d is the 
bead diameter, ryo a dimensionless number and p pc 
has been assumed), so that the local shear rate obeys 
T^/'^ — rioP/{Pc(P'^xy)i- This local relationship can be 
compared to the experimental data by plotting the local 
RMS velocity as a function of the local velocity gradient 
7. As shown in Fig. |^, a remarkable scaling relationship 
is found between these two local quantities, 

yl/2 
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RMS 
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(3) 



with an exponent a ~ 0.4. This scaling behaviour is not 
in agreement with the simple scaling of the viscosity with 
the collision frequency only, for which an exponent of 1 
would have been measured. 

The actual scaling can be understood if we assume 
that the dimensionless coefficient 770 contains an addi- 
tional contribution that diverges algebraically, when the 
density reaches random-close packing, pc, i.e. 770 = 
770(1 — pI Pc)~^ ■ The density can be related to the granu- 
lar temperature through the constant isotropic pressure 
P, which can be formally written as P = pf{p)T. The 
dimensionless function f{p) is expected to diverge like 
(1 — p/pc)^^ at random close packing so that at 
constant pressure P, in a region where p pc, one has 
Pc — P ^ T . Note that this relation is consistent with 
Reynold's remark that the density has to decrease in 
order to allow flow, which occurs for non-zero granu- 
lar temperature. Together with the constant shear rate 
condition, these relations imply a nonlinear algebraic 
scaling of 6V{y) with the shear rate, with an exponent 
a = (2/3+ With the experimentally determined ex- 

ponent a = 0.4, we obtain the overall divergence of the 
viscosity with density as 77 ~ [pc — p)^'^^. 

Using the previously obtained temperature profile, Eq. 
(^) can be integrated (with no-slip boundary conditions 
at both walls) to give the velocity profile. The corre- 
sponding result is plotted in Fig. |^ as a dashed line, 
indicating good agreement with the experimental profile. 
The calculation also yields a relation between the bound- 
ary temperature Tq and the shear rate U, as Tq ~ U" in 
agreement with the experimental data (not shown). 
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FIG. 4. Connection between the local RMS velocity fluc- 
tuations and local shear rate (same symbols as for Fig. 
Local fluctuations are found to increase approximatively as a 
power law of the local velocity gradient, with a power of 0.4 
(dashed line). 

Finally, we measure granular flow when the gap be- 
tween the two cylinders is only 5 particle diameters wide. 
In this case, the inner cylinder jams and can only be 
moved at high airflow rates and high rotation rates. For 
the narrow gap we find that V{y) decreases linearly with 
y and that 6V(y) is uniform across the gap at shear rates 
of 0.1 rev. /sec and 0.02 rev. /sec. This result is consistent 
with our model: A linear velocity profile is expected for 
a constant granular temperature. A constant granular 
temperature in turn is expected for cell widths smaller 
than the previously defined decay length 5. 

The success of the present local hydrodynamic model in 
describing the experimental data (even for small driving 
velocities) calls for a deeper understanding. The basic 
ingredient of the present model is the algebraicly diverg- 
ing viscosity close to random close packing. In a granular 
material, the ability to flow relies mainly on a decrease in 
the local density. Such a decrease can be caused by the 
onset of RMS velocity fluctuations. Lowering the den- 
sity allows some particles to move, creating "fluidized 
regions", while other regions remain in clusters at rest. 
The flow properties can be thus formulated in terms of 
the local fluidization probability, or cluster lifetime. In 
a granular material, these quantities might be estimated 
by means of a free volume estimate, along the lines pro- 
posed by Edwards |12| fo r granular matter and earlier by 
Cohen and TurnbuU |13| in the context of transport in 
glassy systems. These models predict a very strong di- 
vergence of the viscosity close to random close packing, 
in qualitative agreement with the present results. More 
precisely, our data are consistent with an algrebraic di- 
vergence, which is also expected in supercooled liquids 
when the density moves off the critical density [Q. 

The model quantitatively describes all aspects of the 



experimental results: (i) The RMS velocity fluctuations 
increase with the local shear rate to a given power, (ii) 
The velocity decreases strongly away from the moving 
wall, (iii) The normalized velocity profile is independent 
of shear rate and shear dynamics. 

More generally, our experimental results indicate that 
there may be a useful analogy between the dynamics of 
granular materials and the behavior of supercooled liq- 
uids close to the glass transition. Such an analogy was 
proposed in a recent paper by Liu and Nagel |15|. Our 



results support this conjecture and indicate the quan- 
titative relationship: The flow properties reported here 
are quantitatively predicted from a locally Newtonian, 
continuum model, provided that the local temperature is 
identified with (SV)^ and that the viscosity is assumed 
to diverge as the density approaches a critical value. 
The specific properties of individual particles only en- 
ter through the few adjustable parameters of the model. 
Once determined in a specific geometry, the flow proper- 
ties in any geometry can be predicted. Further experi- 
ments are underway to verify this predictive power of the 
model. 
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